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We present a unitary dispersive model for the r; —> Stt decay process based upon the Khuri- 
Treiman equations which are solved by means of the Pasqnier inversion method. The description 
of the hadronic final-state interactions for the ?7 —>■ Stt decay is essential to reproduce the available 
data and to understand the existing discrepancies between Dalitz plot parameters from experiment 
and chiral perturbation theory. Our approach incorporates substraction constants that are fixed by 
fitting the recent high-statistics WASA-at-COSY data for r) —>■ Based on the parameters 

obtained we predict the slope parameter for the nentral channel to be a = —0.022 ±0.004. Through 
matching to next-to-leading-order chiral perturbation theory we estimate the quark mass double 
ratio to be <3 = 21.4 ± 0.4. 
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I. INTRODUCTION 

Production of three particles plays an important role 
in hadron physics. It sheds light on the reaction dynam¬ 
ics, e.g. the OZI rule, and can amplify production of 
hadron resonances, with the mysterious XYZ states seen 
in the spectrum of charmonia and bottomonia [T] being 
the most recent examples. The need for precision analysis 
of final states containing three light hadrons has become 
even more pressing given the high quality data emerging 
from the various hadron facilities around the world, in¬ 
cluding Jefferson Lab, COMPASS and BESIII [5H5]. Re¬ 
cently, significant progress has been made in analysis of 
hadron-hadron interactions at low energies based on the 
S-matrix principles of unitarity, analyticity and crossing 
symmetry [MS- At low energies, unitarity is an impor¬ 
tant constraint given that there is only a limited number 
of contributing channels. Unitarity also determines the 
analytical properties of partial waves and constraints res¬ 
onant scattering. Implementation of crossing-symmetry 
is much more difficult since it is related to the underlying 
dynamics. However, at low energies it can be systemati¬ 
cally investigated by identifying the most important, i.e. 
closest to the physical region, singularities of the cross¬ 
channel amplitudes, and for example in reactions involv¬ 
ing Goldstone bosons these can be constrained by chiral 
symmetry of QGD [lUlfTT] . 

In this paper we focus on decays of the rj meson to three 
pions. From the experimental side, the high-quality data 
from WASA-at-GOSY |I2L I13j . Crystal Barrel [Hill], 
and KLOE [ini E] , along with the data from CLAS [3| , 
which is currently being analyzed, present an opportu- 
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nity for precision analysis of the Dalitz distribution. In 
the charged decay channel, rj —>■ we only have 

access to the binned data from the WASA-at-COSY [12] 
experiment and therefore it is the only data set we use 
in our data-driven analysis. From the theoretical point 
of view ?7 —>■ Stt decays are of interest because of isospin 
violation. These decays are dominated by the intrinsic 
isospin breaking effects in QGD as electromagnetic ef¬ 
fects are expected to be small [Mils]- Consequently, the 
decay width for ry —>■ Stt is expected to be proportional to 
the light quark mass difference and the decay amplitude 
is often expressed in terms of the quantity, l/Q^ defined 
by 

1 ml- ml 
ml — mf 

Here to = (rriu + my) 12 is the average of the u and d 
quark masses. One determines Q by comparing a the¬ 
oretical prediction with the experimental decay width 
r{r] —i 7r+7r“7r°) = 281 ± 28 eV [1]. However, it is im¬ 
portant to emphasize that this procedure requires that 
the amplitude implements chiral constraints or at least 
it agrees with the leading-order chiral perturbation the¬ 
ory (xPT), which is where Q originates. Once Q is ex¬ 
tracted, it can be combined with the knowledge of the to 
and TOs, e.g. from lattice simulations, to determine the 
light-quark mass difference. 

It is necessary to consider the 77 —?► Stt decay ampli¬ 
tudes beyond yPT. This is apparent when considering 
contributions to r (77 —> tt+tt^tt*^) from the first few terms 
in the low energy expansion. Specifically, the leading- 
order yPT result, = 66 eV [SU] El], is ap¬ 

proximately four times smaller than expected. Inclu¬ 
sion of next-to-leading (one loop) corrections increases 
the theoretical prediction to = 167 ± 50 eV 
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\n \ , which is still significantly below the data. The next- 
to-next-to-leading calculation (two loops) has been per¬ 
formed recently [23] ■ It pushes the decay width further 
towards the data; however, it contains a large number 
of low energy constants. In addition to the apparent 
poor convergence, low orders of xPT give an incorrect 
result for the shape of the Dalitz distribution in the 
neutral 37r° decay. To the leading order, this distribu¬ 
tion is represented by a single parameter, a, which xPT 
predicts to be positive while the experimental result is 
a =—0.0317 ± 0.0016 [T]. The fact that chiral expan¬ 
sion converges slowly indicates the importance of final 
state interactions. This is expected to be a consequence 
of unitarity, which in xPT is incorporated only order by 
order. To fulfill unitarity various dispersive frameworks 
were developed PH [23] with recent updates of [21123 
and [21] . These analyses are based on the Khuri-Treiman 
(KT) representation [29]. In the KT approach, partial 
waves are given in the elastic approximation with the left- 
hand cut contributions computed from cross-channel am¬ 
plitudes that are approximated by the same elastic par¬ 
tial waves as in the direct channel and are bootstrapped. 
Other calculations employed, for example, nonrelativis- 
tic effective field theory (NREFT) [30] and alternative 
dispersive approaches were studied in [3T] . 

The final state interactions in 77 —>■ 37 r at low en¬ 
ergies can be approximated by elastic tttt scattering. 
These amplitudes are available with high precision up 
to ^/s = 1.1 GeV [ 3 . However, dispersion calculations 
involve integrals over all energies. In order to suppress 
the unknown high-energy region, the dispersive integrals 
are usually over-subtracted and the subtraction constants 

are fixed by comparing to the data [23 [32]- In [33] the 

authors used an alternative method whereby the disper¬ 
sive integral was split into elastic and inelastic contribu¬ 
tions and the latter was parametrized by a power series 
in a suitably chosen conformal variable. In the current 
work, we apply yet a different approach. We obtain the 
solution of the KT equation using the so-called Pasquier 
inversion method [331133] • In this case the dependence on 
the unknown high energy region is traded for by the de¬ 
pendence on the far left-hand cuts. The advantages and 
disadvantages of alternative procedures were discussed in 

m- 

The paper is organized as follows. In Section [^ we 
present the basic formalism and discuss how three body 
effects are incorporated using the Pasquier inversion. The 
numerical results are presented in Section m which we 
divide into two parts. In the first part we perform a data- 
driven dispersive analysis of the WASA-at-COSY data 
m without input from xPT. We show the fitted Dalitz 
plot parameters for the charged decay and predict the 
slope parameter for the neutral decay channel. In the 
second part we match our amplitudes to xPT in order 
to extract the Q value. Conclusions are summarized in 
Section |IV[ 


II. FORMALISM 

A. Kinematics and partial wave expansion 

The isospin violating 77 —> 37 r decay involves a 

AJ = 1 interaction. The transition matrix elements, 
t, m), depends on four isospin indices, with the 
index 77 referring to the isospin component of the in¬ 
teraction and a, /3 ,7 to three pions. In terms of the 
particle momenta the three Mandelstam variables are 
s = {pi +P 2 f = iPi-Psf, t = {P2+P3y = (P4 -Pl)^ 
and u = {pi + pa)'^ = (P 4 — ^ 2 )^- The Mandelstam vari¬ 
ables satisfy s -|- t -I- u = + ttt,^ + m\ -\- in\, with tti^ 

being the mass of the 77 , also referred to as particle i = 4 
and rrii, i = I..3 to the pions. On account of crossing 
symmetry, the following processes are described by the 
same complex function (with the bar denoting an an¬ 
tiparticle): the s-channel scattering 4 -f 3 —)■ I -I- 2, the 
t-channel scattering 4 -f 1 —)■ 2 -|- 3, the u-channel scatter¬ 
ing, 4 -f 2 —)■ 1 -f 3, and the decay channel A —>■ 1 + 2 + 3. 
In particular the amplitude in the decay channel will 
be derived by analytical continuation of the s-channel 
partial wave expansion. In the s-channel, the amplitude 
A“^^’’(s, t, m) has the following partial wave (p.w.) de¬ 
composition, 

00 

A“^^''(s,t,77) = ^ ^(2L + l)PL{zs)vi^j^^AjUs), 

L^O I 

( 2 ) 

where Pl{zs) is the Legendre polynomial and is a co¬ 
sine of the center-of-mass scattering angle Og, 

s{t-u) + {ml - ml) (m^ - mj) 

8 COS s A1/2 (s, 777^, 77r|) A1 A(s, 777.^, 77r|) 

The usual Kallen triangle function is given by 
X{a,b,c) = + b'^ + — 2 {ab + b c + a c) and (/, L) la¬ 

bel isospin and orbital angular momentum quantum 
numbers in the s-channel with I + L = even due to Bose 
symmetry of pions. The isospin projection operators 
T’ap'-ir] given in Appendix A We note that at this 
stage the partial waves are arferarily normalized. The 
unitary relation, which we discuss in the following is ho¬ 
mogeneous in A and at the end we will normalize the 
amplitude by comparing with the experimental data. 

The p.w. amplitudes A/i(s) have both the right-hand 
cut discontinuities demanded by the direct channel uni¬ 
tarity and left-hand cut discontinuities from exchanges 
in the t and u channels. We emphasize that Eq. ([^ is 
exact in the s-channel physical region, when the infinite 
sum over L converges. The amplitudes in the other chan¬ 
nels are obtained by analytical continuation. Low-energy 
approaches based on partial wave expansion involve trun¬ 
cation of the partial waves series at some L — Lmax < 00 , 
which violates analytical properties of cross-channel am¬ 
plitudes. To partially recover those, we represent the 
amplitude as a sum of truncated partial wave series in 



3 



each of the three channels [laiMlISTHlI], 


5 ]( 2 L + 1 ) 

L^O I 

+PL{zt) Vp^ar,aiL{t) + Pl{zu) l{u)^ , 

(4) 


where the amplitudes are qjl defined as having only 
right-hand discontinuities demanded by unitarity in the 
respective channels. The center of mass scattering angles 
in the t- and the rt-channel are given by 

t{s — u) + (to§ — m|) (to^ — m\) 

TO 2 , m§) 

_ M (t - s) -f {ml - mf) {ml - m|) 

We remark that the decomposition in Eq. 0 satisfies 
crossing symmetry explicitly; however, violation of an- 
alyticity remains since the amplitude contains a finite 
number of high-spin partial waves in any given channel. 
This would be a problem at high energies but hopefully 
does not influence our low-energy analysis. What the 
representation in Eq. Q does is to allow for unitarity 
to be implemented in all three channels. We also note 
that decomposition in Eq. 0 is exact up to NNLO in 
yPT [m 03] and is often referred to as “reconstruction 
theorem”. 

It is convenient to express the p.w. amplitude Ajl(s) 
{c.f. Eq. 0 ) in terms of the amplitudes 0 /^( 5 ) that are 
defined by Eq. Q , 

Aj,is) = afl^^Hs) + a^l^\s). ( 6 ) 

Here the amplitude has only the right-hand dis¬ 

continuity, 

=aiLis), (7) 

and the left-hand discontinuities of originate 

from the exchange terms, 


^Left 

^IL 


W=EE 


{2L' + 1) 


dZsPL{Zs) 


L'=o r 

'll!' 


/-I 


X {PL'{zt)Cll'avL'it) -k PL'{zu)Cl^^arL'{u)). 


( 8 ) 


B. Unitarity and the three-body effects in the 
decay channel 


In the following we consider both decay modes of the 
r] meson, the charged decay 77 —>■ 7 r+ 7 r“ 7 r°, and the neu¬ 
tral decay 77 —>■ When comparing with experimen¬ 

tal data it is important to have an accurate descrip¬ 
tion of the phase space boundary, thus in the compu¬ 
tation of the kinematical factors we use the physical pion 
masses. Elsewhere we assume the isospin limit and use 
rrii = (2m,n.+ -I- TO7ro)/3 = i.e. the isospin averaged 
mass. 

The model is defined by Eq. ([^ together with the 
elastic unitarity constraint for the right-hand disconti¬ 
nuity mi, 

^ Right, ^ ^ f Right, ■ , _ Right, _ - A 
^ajL aj^ {b it) j 

= fids) p{s) {afd'^^s) + ddh^s)), (9) 

where p{s) = a/ 1 — dm^/s. The elastic tttt partial 
wave amplitudes are denoted by //l and normalized by 
Im(l////,(s)) = —p{s). Therefore, the amplitudes aiL{s) 
satisfy the relation. 


^aids) = fkis)p{s)(ajL{s) + Y, 


[ dtPL{zs)PL'{zt)Cdt arL'{t)\- 

Jt_{s) J 


( 10 ) 


The first term on the right-hand side of Eq. (10) rep¬ 


resents the contribution from the direct s-channel, 
4 -I- 3 —>■ 1 -I- 2, to the s-channel partial-wave projection 
of the unitarity relation and it is illustrated in the di¬ 
agram in Fig. [^a). The second term, illustrated in 
Fig.^b), gives the contribution from the exchange con¬ 
tributions in the t-channel 4 -|- 1 —> 2 -|- 3, and rt-channel 
4 -I- 2 —>• 1 -I- 3. In Eq. (101, using Eq. ([^, we changed the 
integration over the Zs to integration over t. 


dzc 


(...) = 


dt 


P+U) 

Lis) K{s) 


S {■■■): 


( 11 ) 


with the integration limits t±(s) corresponding to 

-Zs = ± 1 , 


Here Cgt and Cs„ are the standard crossing matrices and 
are given in Appendix [X] 


t±{s) = 


m, 


^ -I- 3 rrii — s 


± 


K{s) 


( 12 ) 


2 


2 s 
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t^(s) 



FIG. 2: Integration contour in the complex t plane. The 
arrows indicate the direction of increasing s in the interval 
from 4 to oo. The points labeled a through i correspond to 
specihc values of s, with (a) t_(oo) = 0, (b) + = 

m^(m^ - rur,), (c) = m^(m^ + m^), (d) 


t-C 


'-) = 4m^, (e) t±(4m^) = 


(f) + 


m^)) = (m^ - m^) , (g) t+((m^ - m^) ) = m^(m^ + m^), 
(h) t+{{mn + m^)^) = m^{m^ — m^), and (i) t+(oo) = —oo, 
respectively. 


The Kacser function K{s) is given by the product of the 
triangle functions and has the following determination 

[MSS] 


{ +k(s) , 4 m^ < s < (m^ — 

i k{s) , (to,, — rriTr)'^ < s < (to,, + to„)^, 

— k(s) , (to,, + TO,r)^ < S < + 00 , 

k(s) = |A(s, ml, ml) X{s, to^, to^)| i /2 . ( 13 ) 

In the scattering region s>(TO,, + TO,r)^ the inte¬ 
gral in Eq. 0 is well defined; however, when 
4 to^ < s < (to,, -l-TO.n')^, analytical continuation to the 
decay region is needed. For this a positive infinitesimal 
imaginary part is added to the eta mass IMBSIESI, 
which leads to the integration contour in the t-plane 
shown in Fig. It is worth noting that the contour 
avoids the unitary cut. Finally, the amplitudes 0 /^( 5 ) 
are obtained by bootstrapping the dispersion reaction 


1 

aiL{s) = - ds 
J 4m2 


, Aa/L(s') 


S' — s 


(14) 


with aiL appearing on the right-hand side {cf. Eq. 0 ) 
together with the input two-body scattering amplitudes, 
fiL{s). 

As in the standard N/D approach, the inhomogeneous 
part in Eq. ( [l0| can be accounted for writing a/L(s) as 
a product of //l(s) times another function of s, whose 
discontinuity is given by the s-channel projection of the 
cross-channel amplitudes. It is also convenient to remove 
any zeros of //l(s), e.g. the Adler zero, since these are 
process dependent. Einally, the partial waves have kine- 
matical singularities, which do not contribute to the dis¬ 
continuity relation given by Eq. (10). Thus, we write 


(15) 


where the first factor removes the kinematical singulari¬ 
ties 


Zl{s) = 


Kis) 


n L 


s/4: — 


(16) 


and the second factor removes zeros from the tttt ampli¬ 
tude, 

, (s-4^V(s-Sa^)) L = 0 , 

^il{s) = { (17) 

1, L>0. 

That is, we assume fiL has zeros in the S'-wave only. 
Note that at leading order in yPT, Adler zeros are lo¬ 
cated at = ml /2 and = 2 ml in the tttt S-wave 

isoscalar and isotensor amplitudes, respectively, and at 
= 4/3 TO^ for r] Stt. In the actual calculation we 
use as input the tttt amplitudes from the phenomenolog¬ 
ical analysis of [7] which have zeros at the same position 
as the leading order in yPT; when matching 77 —> Stt 
with yPT we use NLO calculation which places the ze¬ 
ros in 77 —Stt at = 1.25 to^ and = 2.7 ml in the 
isoscalar and isotensor channels, respectively. 


Finally, it follows from Eq. (10) and Eq. (15) that the 


function giL das the discontinuity given by 

Ac 7 /l(s) = - 6 »(-s) giLis) 

fids) 

+ e{s - Ami) E Pis)PL{zs) 

L '=0 /' 


(18) 


K{s)/s Fil{s)Zl{s) 


ft+is) 

X / dt PL'{zt)Cd ZL'{t)J^i'L'{t)fl'L'{t) gi'L'{t) . 


The first term on the left-hand side takes into account the 
left-hand cut of fids)] i-e. in addition to the unitary cut, 
giL has a left-hand cut determined by fiL to guarantee 
that there is no dynamical left-hand cut in the amplitudes 
ajL- The integrand in Eq. (18) is free from kinematical 


singularities in t and the function giL{s) satisfies 

, /XgiL{s') 


1 p Ag, 

gids) = ds —j 

TT 7-00 S 


(19) 


Inserting Eq. (18) into Eq. (19) we obtain a double in¬ 


tegral equations for giL{s), which can be reduced to a 
single integral equation by changing the order of disper¬ 
sive integral (over s) and the angular projection (internal 
over t). The procedure, which we referred to earlier as the 
Pasquier inversion, was developed in jMISSj and recently 
revisited in [^. It leads to the following representation 


gids) = — 


/■« ,, 1 AfiUs) 


1 n{M — mTr)^ Lmax 

~ E ^^iLJ'L'{s,t) 

L '=0 /' 


X Cli fl'L'if) gi'L'd), (20) 


aiL^s) = Zl{s)Pil{s) fiL{s) giL{s) 
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where the kernel function (s, t) is given explicitly 

in Appendix [ b| The left-hand cut contribution to giL^s) 
is largely unknown. Since we are primarily interested in 
the physical decay region we therefore parametrize con¬ 
tributions to giL from integration over s < 0. In the 
simplest approximation these are reduced to a constant. 
A more elaborated representation could, for example, in¬ 
volve a conformal map of the s-plane cut along the neg¬ 
ative real axis onto a unit circle m- However, in the 
analysis of the data we find the simple approximation to 
be sufficient: 


1 /*(M — Ljria x 

gids) = gihiso) + - 

L'=0 I' 

X {K.il,I'L' (s, t) — ICjlJ'L' {sq, (t) gi'L' (t) ■ 

( 21 ) 


This equation can now be solved using standard matrix 
inversion methods with the subtraction constants giL^so) 
as fitting parameters. The subtraction point is arbitrary 
and we choose it to coincide with the Adler zero of the 
LO xPT So = 4/3 . After solving the integral equation 

for giL{s), we compute aids) from Eq. Finally, 

to compare with the experimental data we convert the 
isospin amplitudes to the charge amplitude, A^{s,t,u) 
for the r] —>■ 7r+7r“7r° and A^(s, t, u) for the neutral case. 
These are given by Eq. Q, 


A^{sA,u) = ^ 


{2L + 1) 


L=0 


g Pl{zs) (aoLis) - a2L(s)) 


+ PLdt) {aiLid +a 2 L{t)) - Pdzu) (aiiiu) - a 2 L{u)) 


A^{s,t,u) = ^ 


{2L + 1) 


L=0 


Pl{zs) (aoL(s) + 2a2L(s)) 


-f (s —^ t) -h (s —y u) 


( 22 ) 


III. NUMERICAL RESULTS 


a.k.a. final-state interactions in the decay region, by per¬ 
forming two analyses. In the first, we do not include 
cross-channel effects and approximate gid^) fo Eq. (21) 
by a constant, setting gid^) = gid^o)- ^ corresponds 
to a traditional isobar model, but with a fully incorpo¬ 
rated two-pion interaction. In the second, we include 
cross-channel rescattering effects and solve Eq. (21). In 
the following we refer to the two cases as “two-body” and 
“three-body”, respectively. 


A. Fitting WASA-at-COSY data 

1 . ?7 —>■ 7r^7r“7r° 

In this subsection we summarize the results of the fit 
to the recent WASA-at-COSY data on ry —>• 7r+7r“7r° [T^ . 
where binned Dalitz plot is given. Up to a normalization 
factor, the Dalitz plot distribution is given by the ampli¬ 
tude squared. 


(23) 

It is convenient to express the amplitude in terms of two 
independent, dimensionless variables {x, y) which are lin¬ 
early related to the Mandelstam variables by 


V3 


y = 


2 ttZyi c 

3 

2 ttZyi c 


{t-u), 

{{rrir, - m^o)^ - s) - 1, 


(24) 


where Qc = to,, — 2 to+ — to° (for the neutral decay we 
use Qn = rrij^ — 3to°). A general property of these vari¬ 
ables is that the physical region of the Dalitz plot lies 
inside the unit circle -I- < 1 centered at a; = y = 0. 

We fit our model to the data |12] by minimizing the 
defined by 


N 

E 

bins 



-|A^({y,^(so)})P 


2 


(25) 


In this section we present our results for the decays 
r] —>• 7r+7r“7r° and y —37r°. We study the systematic un¬ 
certainties of the model by using different sets of partial 
waves, i.e. varying Lmax and maximal isospin. We have 
found that partial waves with (L > 2) are negligible in 
the physical decay region, 4 to^ < s < (to,, — to,,)^. As 
input we use two-pion scattering amplitudes from the 
analysis of [7]. The parameters of the fit are the sub¬ 
traction constants, y/L(so), for each contributing partial 
wave. Our aim is to fix these by fitting y —>■ 7r"''7r“7r° 
decay using the high statistic WASA-at-COSY data [T^ 
and by matching to NLO yPT [22]. The results for the 
y —>■ 3d decay mode will then constitute a prediction, 
which we compare with the Dalitz plot distribution from 
[48| . We investigate the role of cross-channel exchanges. 


over the set of subtraction constants, gid^o)- fo 
Eq. (25), is the acceptance-corrected number of 

events m each of the N = 59, Ax = Ay = 0.2 wide mass 
bins. The data is normalized to unity at x = y = 0 
and A|A Lata is the statistical uncertainty. Note, that 
since Eq. ( |21| ) is linear in y/^, the parameter 500 ( 50 ) can 
be factored out and fixed by the overall normalization. 
Since normalization of the data is arbitrary the absolute 
value of 500 ( 50 ) is irrelevant. Therefore, in Tablejij which 
summarized fit results, when presenting results of two- 
body fits we quote {gjido) ± ^gjldo))/god^o)■ When 
presenting results of three-body fit we quote ( 57 ^( 50 ) ± 
^ff/L(5o))/3oo(so), where 5oo(5o) is the central value ob¬ 
tained in the two-body fit with the same number of par¬ 
tial waves. We do the latter to illustrate the relative 
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TABLE I: Results of two-body and three-body fits for different wave sets. 



ffoo(so)/5oo*’^ 

ff2o(so)/5of^ 

S'ii(so)/Soo'’’ 

X^ld.o.f. 

(7,L) = (0,0) 

two-body 

1.000 ±0.002 

- 

- 

2.2 

three-body 

1.062 ±0.002 

- 

- 

15 

(/,L) = (0,0), (2,0) 

two-body 

1.000 ±0.003 

0.04 ±0.01 

- 

1.69 

three-body 

1.138 ±0.003 

0.29 ±0.01 

- 

1.67 

(7,L) = (0,0), (1,1) 

two-body 

1.000 ±0.002 

” 

0.058 ± 0.009 

1.45 

three-body 

1.043 ±0.005 


0.233 ± 0.009 

0.95 (Set 1) 

(7,L) = (0,0), (2,0), (1,1) 

two-body 

1.00 ±0.02 

-0.26 ±0.05 

0.38 ±0.07 

0.94 

three-body 

1.19 ±0.01 

0.14 ±0.03 

0.28 ±0.04 

0.90 (Set 2) 

TABLE 11: Dalitz plot parameters for y —>■ 7r"'"7r 7r°. 
(7,L) = (0,0), (2,0), (1,1) cases respectively (see Table|^. 

Set 1 and Set 2 correspond to (7, L) = 

(0,0), (1,1) and 


a b 

d 

/ 
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WASA-at-COSY [H] 

-1.144 ±0.018 0.219 ±0.019 ±0.037 0.086 ± 0.018 ± 0.018 0.115 ±0.037 

- 

KLOE [16] 

-1.090 ±0.005lH?^ 0.124 ± 0.006 ±0.010 0.057 ± O.OOOt^J));^ 

0.14 ±0.01 ±0.02 - 

CBarrel [T3] 

-1.22 ±0.07 0.22 ±0.11 

0.06 ± 0.04 (fixed) 

- 

- 

Layter et al. [49] 

-1.080 ±0.014 0.03 ±0.03 

0.05 ±0.03 

- 

- 

Gormley et al. |5U| 

-1.17 ±0.02 0.21 ±0.03 

0.06 ±0.04 

- 

- 

Theory 

Set 1 

-1.116 ±0.030 0.188 ±0.010 

0.047 ±0.005 

0.093 ± 0.004 

-0.020 ± 0.006 

Set 2 

-1.117 ±0.035 0.188 ±0.014 

0.079 ± 0.003 

0.090 ± 0.003 

-0.063 ±0.012 

NLO [22] 

-1.371 0.452 

0.053 

0.027 

- 

NNLO [23] 

-1.271 ±0.075 0.394 ±0.102 

0.055 ±0.057 

0.025 ±0.160 

- 

Kambor et al. [24] 

-1.16 0.24...0.26 

0.09...0.10 

- 

- 

NREFT [3D] 

-1.213 ±0.014 0.308 ±0.023 

0.050 ± 0.003 

0.083 ±0.019 

-0.039 ± 0.002 


change in normalization between two- and three-body 
fits. 

In the first fit we use a single, scalar-isoscalar, ago par¬ 
tial wave. In this case, the model gives a parameter free 
prediction for the event distribution. We observe that the 
(I, L) = (0, 0) amplitude provides the dominant contri¬ 


bution that covers approximately 90% of the Dalitz plot. 
The calculated x^/d.o./. for the two-body and three- 
body cases are 2.2 and 15, respectively. In Fig. (upper 
panels) we compare our results and the data projected 
onto the x and y axes. The error bars associated with 
the model originate from the uncertainties in the pion- 
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pion amplitude fiL [Z] and from the statistical error in 
fitting the overall normalization. 

In the next step, we add the isospin-2 S'-wave. In this 
case we fit two parameters, one gives the overall normal¬ 
ization and the other contributes to a modification of the 
shape of the Dalitz plot. The resulting parameters and 
X^/d.o.f are given in Table]^ In both, the two- and three- 
body fits we find that the model slightly underestimates 
the data. The inclusion of the second (/, L) = (2,0) wave 
significantly improves and also drastically reduces the 
difference in the fit quality between the two- and three- 
body cases pertinent in the fit with the single {I, L) = 0 
wave. 

In the spirit of keeping the number of free parame¬ 
ters as low as possible, we considered another set of two 
waves, (/, L) = (0, 0), (1,1), before taking into account a 
complete sum of S and P waves. In this case there is also 
one parameter that affects the shape of the Dalitz distri¬ 
bution and we find x^/d-o.f = 1.45 and x^/d-o.f = 0.95 
in the two-body and three-body fits, respectively. Hence, 
it seems that the data favor the isovector P-wave contri¬ 
bution over the isospin-2 S'-wave. The results of the fit 
are shown in Fig. 

We now turn to the case when a complete set of S and 
P waves is incorporated, i.e. (I, L) = (0,0), (2, 0), (1,1). 
The two- and three-body fits result in a comparable 
x'^/d.o.f around 0.9. 

It is instructive to compare the results of the three- 
body fits. In the fit with a single (/, L) = (0,0) am¬ 
plitude, the three-body fit converges poorly indicating 
importance of higher partial waves that are brought in 
by the cross-channel exchanges. Thus apparent conver¬ 
gence of the two-body fit in this case is deceptive. With 
any combination of higher partial waves all calculated 
three-body x^/d.o.f are quite similar to the two-body 
fits, except for the case when only (/,T) = ( 0 , 0 ), ( 1 , 1 ) 
amplitudes were considered. 

Often, an effective range expansion of the Dalitz plot 
near x = y = 0 is used to parametrize the y decay distri¬ 
bution. For the charged decay it leads to 


|HC(0,0)P 


1 + ay + by^ + cx + dx^ 

+ exy + fy^ + gx'^y -. 


(26) 


The charge conjugation symmetry, x —>■ —x requires 
terms odd in x to vanish, i.e. c = e = 0. In Table 
[Tl| we give the Dalitz plot parameters from our three- 
body fits based on the (/, T) = ( 0 , 0 ), ( 1 , 1 ) (set 1 ) and 
(/, L) = (0,0), (2,0), (1,1) (set 2) wave sets. For com¬ 
parison we quote the results of next-to-leading-order 
(NLO) and next-to-next-to leading order (NNLO) of 
yPT [53] , the dispersive analysis from [53| , NREFT 

[30j and alternative dispersive approach m- We also in¬ 
clude Dalitz parameters extracted from direct fits to the 
experimental data na nn [Miiisiiio]. The most recent 
analyses where performed by the WASA-at-COSY [T5] 
and KLOE [T3] collaborations. As expected, our Dalitz 


plot parameters are consistent with the WASA-at-COSY 
parameters within the error bars. We also observe that 
central values of the fit tend toward the KLOE results. 


2. rj ^ 37 r° 


The results obtained in the charged mode can be 
used to predict the Dalitz plot parameters for the neu¬ 
tral channel. The Dalitz parameters are defined as 
coefficients in the expansion around the center of the 
Dalitz plot using the polar coordinates x = y/z cos 4> and 
y = y/z sin ([> in Eq. (24) 


|A^(0,0)P 


= 1 -\-2az + 2 (3 sin + ■ ■ 


(27) 


The slope parameter a has been extracted from several 
experiments, while to the best of our knowledge, there 
is no determination of /3 or higher moments. In Ta¬ 
ble jlllj we compare our findings with the experimental 
measurements and other theoretical predictions. The av¬ 
erage of experimental results compiled by the PDG is 
a = -0.0317±0.0016 [I]. 

As in the case of the charged mode, our results ob¬ 
tained with the the two sets of waves are quite similar. 
The predicted slope parameter is a(Set 1) = —0.023 and 
a(Set 2) = —0.020. Even though both sets describe the 
charged data well, the predicted slope parameter in the 
neutral case is above the PDG value. As shown in [531130] 
the Dalitz plot parameters of the neutral and charged de¬ 
cays are related by 


a = (d + b-\a^- Im(a)^ 

4Q2 4 


" 4Q2 


< 


d + b — \ 
4 


(28) 


where the factors Qc, Qn were defined below Eq. (24|. 
Note that we only take Qc Qn in the overall normal¬ 
ization while we use Qc = Qn when solving dispersion 
relations for the partial wave amplitudes. Here, the com¬ 
plex parameters a is the coefficient of the linear term in 
the expansion of the charged amplitude A'^{x,y), 


{x, y) oc 1 + ay + ... 


(29) 


Using the Dalitz plot parameters from WASA-at-COSY 
and KLOE collaborations one finds 


,WASA < _o.oo6 , < -0.033. (30) 


The large difference in the upper limits is due to the 
difference in the b parameter which differs by a factor 
of two between the two data sets. As pointed out in 
[30] the value for Im(a) can be sizable due to tttt fi¬ 
nal state interactions. Our results confirm this find¬ 
ing and we obtain Im(a) = —0.18 ±0.03. Nevertheless, 
since = —0.006 is quite large the Im(a) term 









FIG. 3: Upper and middle panels are the x- and y-projection plots. Black circles are the data. Red squares and blue squares 
represent results of the two-body and three-body fits, respectively. The hts are performed on the Dalitz distribution m sown 
in the bottom left panel using a single, (7, L) = (0,0) wave (upper panels) and two waves, (7, L) = (0,0), (1,1) (central panels). 
For better visualization fit results are shifted horizontally (three-body to right and two-body to left) from the experimental 
points. The bottom right panel is the Dalitz distribution from the three-body fit with (I,L) = (0,0), (1,1) waves. 
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TABLE III: Dalitz plot parameters for rj —>■ Srr®. Set 

1 and Set 2 correspond to (7, L) = (0,0), (1,1) and 
(7,7/) = (0,0), (2,0), (1,1) cases respectively (see Table]^. 


a 


GAMS-2000 [H] 

-0.022 ±0.023 

- 

Crystal Barrel, LEAR |52| 

-0.052 ±0.020 

- 

Crystal Ball, BNL [IS] 

-0.031 ±0.004 

- 

SND [53j 

-0.010 ±0.023 

- 

CELSIUS-WASA [B] 

-0.026 ±0.014 

- 

WASA-at-COSY |5l| 

-0.027 ±0.009 

- 

MAMI-B [53] 

-0.032 ± 0.004 

- 

MAMI-C [IS] 

-0.032 ± 0.003 

- 

KLOE [T7] 

-0.0301 ±0.0050 

- 

PDG average [1] 

-0.0317 ±0.0016 

- 

Theory 

Set 1 

-0.023 ±0.004 - 

-0.000 ±0.002 

Set 2 

-0.020 ±0.004 - 

-0.001 ±0.003 

NLO [H] 

±0.013 

- 

NNLO [23[ 

±0.013 ±0.032 

- 

Kambor et al. [23] 

-0.007...-0.014 

- 

NREFT [30] 

-0.025 ±0.005 - 

-0.004 ±0.001 

Kampf et al. |31| 

-0.044 ± 0.004 

- 


alone can not be responsible for lowering a to the PDG 
value. Once the KLOE data become available [SSj it 
would be very interesting to perform a combined fit of 
the WASA-at-COSY and KLOE measurements. 

The neutral channel does not depend on the P-wave 
amplitude contributing to the charged decay mode and 
it contains only even partial waves. Unfortunately, using 
the charge mode we could not find sensitivity to the D- 
wave which was omitted from the Table |l] Finally in 
Fig|^ we compare our results with the recent MAMI-C 
measurement [48) . The R{z) function is determined as 

( 31 ) 

fo d(l)9{ip{s,t,u)) 

where 

(p{s,t,u) = stu - mlo - mlo)^ = 0 (32) 

defines the boundary of the Dalitz plot distribution and 
6{x) is the step function. We observe that a cusp around 



FIG. 4: Comparison of R{z) plot from [4^ (black points) with 
our predictions from Table |TTT] that correspond to Set 1 (blue 
band) and Set 2 (red band). 


z ~ 0.765 appears in R{z) for nonzero /3. This is a kine- 
matical effect which reflects the fact that for larger z the 
phase space distribution in the Dalitz plot is no longer 
circular. We find our results for Sets 1 and 2 provide a 
satisfactory agreement with the data. 


B. Matching to xPT and the Q-value 

We remind that the data in |12j were normalized to 
the center of the Dalitz plot and therefore our model 
only predicts the Dalitz plot distributions for the charged 
and neutral decays. The overall normalization can be 
fixed by comparing the experimental decay widths with 
the phase space integral over the corresponding squared 
amplitudes, 


r = 



\A{x,y)\^ 

|A(0,0)P ’ 


(33) 


with the boundaries of the integral determined by the 
phase space. We emphasize that the quantity defined 
in Eq. 0 enters into the normalization constant N. In 
order to determine one has to match the model, dis¬ 
persive amplitude, with xPT where is defined. 

As discussed in Sec. [Tj the xPT [23] series seems to 
converge rather slowly and the question arises to which 
order of the yPT should one match the model. It would 
be desirable to find a matching point where on the yPT 
side contributions, from powers of Mandelstam invari¬ 
ants, are small. Therefore, matching the amplitudes in 
the physical region may not be the best option. Up to 
NNLO the chiral amplitude satisfies the decomposition of 
Eq. and up to this order matching is simplified since 
it is sufficient to match the single variable, partial wave 
amplitudes a/L(s). The xPT amplitude for the charged 
decay, up to NNLO can be written in the form 




1 


M{s,t,u), 


(34) 
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FIG. 5: Upper panels: x- and t/-projections of the Dalitz plots. Black circles represent the data. The red squares and blue squares 
are model results using amplitudes with only two-body and including three-body correlations, respectively. The amplitudes 
were computed using three partial wave components with (I,L) = (0,0), (2,0), (1,1). For better visualization fit results are shifted 
horizontally (three-body to right and two-body to left) from the experimental points. Bottom panels: The comparison of the 
NLO xPT amplitudes Mj’s (black curves), with the two-body (red curves) and three-body (blue curves) dispersive amplitudes. 
Real parts are shown with solid lines and imaginary with dashed lines. In all figures the unknown couplings were fixed by 
matching to NLO xPT (see Eq.(38l). 


TABLE IV: Values of Q from different calculations. 


Theory 

Q 

Set 1 

21.7 ±0.4 

Set 2 

21.1 ±0.4 

Lattice {Nf = 2 + 1)°' 

22.6 ±0.9 

NLO dH 

20.1 

NNLO [231 

22.9 

Kambor et al. 

22.4 ±0.9 

Kampf et al. [31] 

23.1 ±0.7 


“Here and in the following we combined in quadrature the errors 
quoted in EZl- 


where = 92.3 MeV is the pion decay constant and 

M{s,t,u) = Mo(s) - ^M 2 (s) -I- M 2 {t) + M 2 {u) 

-b (s — u)Mi(t) + (s — t)Mi{u). (35) 


Explicit expressions for the functions Mj at various or¬ 
ders in the chiral expansion can be found in |23j . Com¬ 
paring Eq. (1^ and Eqs. (34), (1^ one finds 


«oo(s) = 3 N^pt Mo{s), 

« 2 o(s) = 2 N^pt M2{s), 

aii(s) = ly^xPT —^-^1(5): 
o s 

where 

1 TOl-(m^-TO^) 


(36) 


(37) 


The NNLO yPT calculation was performed in |23) . 
The order 0(p®) LECs were estimated using a reso¬ 
nance saturation model and error analysis was not pro¬ 
vided. Given that uncertainties in the low energy con¬ 
stants entering M/’s at the NNLOs are not quantita¬ 
tively settled in the following we choose to match our 
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dispersive calculation with the NLO xPT result. In 
this case one can use the NLO relations between decay 
constants and meson masses which reduces the number 
of low energy constants in the chiral amplitude to one, 
L 3 = (—2.35 ± 0.37) • 10“^ [SS]. We choose the matching 
point to coincide with the subtraction point in Eq. (21 1 , 
which in turn was chosen to coincide with the Adler zero 
in the LO yPT amplitude. In that case the determined 
parameters from matching are the same for the two-body 
and three-body scenarios. 

In the following we consider two methods for matching 
the dispersive analysis with yPT. In the first case we use 
Eq. (37) together with the xPT NLO amplitudes M/’s 
to compute the overall normalization and the parameters 
57 l(so)) which in turn completely determine dispersive 
amplitudes of our model. We find, 


5oo(so) — 16-1 N^pt , 

goo{so)/92o{so)/9ii{so) = 1/0.12/(0.129 ± 0.014). (38) 


This confirms that the amplitude (/, L) = (0, 0) is domi¬ 
nant. In the lower panel of Fig. we compare the xPT 
amplitudes with the dispersive ones, the latter obtained 
using the subtraction constants from Eq. (38 1 . Compar¬ 
ing with the WASA-at-COSY data shown in the upper 
panel in Fig. j^we find that the dispersive amplitude 
fixed by Eq. ( [^ gives /d-o.f. of approximately 13.0 
using only two-body amplitudes, which is reduced to 2.9 
when three-body rescattering contributions are included. 
Even though the model compares reasonably well with 
the data, the large value of x^/d.o./. prevents us from 
extracting the Q-value using this method. 

To extract the Q-value we therefore use the xPT ampli¬ 
tudes to determine the overall normalization only, while 
for the subtraction constants gipiso) we use the re¬ 
sults from the fit of the WASA-at-COSY data described 
in the previous section. We find (5(Set 1) = 21.7 ± 0.4 
and (5(Set 2) = 21.1 ± 0.4 for the two sets of parameters 
given in Table |TT] Comparison of our findings with pre¬ 
vious results is summarized in Table IIVI We observe that 
the extracted Q-values are somewhat smaller compared 
to [2211H1I3I]) and within Icr from the recent {Nf = 2-|-l) 
lattice computations m- We note that lattice calcu¬ 
lations of electromagnetic correction tor Nf = 2-1-1 
are not yet available, while tor Nf = 2 these were re¬ 
ported in [? ]. The lattice result given in Table IV de¬ 
pends on the input value for the light quark mass ratio, 
iTiu/md = 0.46 ± 0.03 which is the LO xPT result re¬ 
duced by a factor of 8(4)% chosen as an estimate of the 
correction from higher-orders chiral effects m- Alterna¬ 
tively, using the extracted Q-value and the Nf = 2 + 1 
lattice result for ms/rh = 27.46 ± 0.44 m we can esti¬ 
mate mu/md- We find 


— = 0.42 ± 0.02 (39) 

rud 

as an average between Sets 1 and 2. Another useful quan¬ 
tity that can be calculated from our Q and ms/rh is the 


so-called i?-value given by 



= 32.2 ±1.3. (40) 

IV. CONCLUSIONS 

In this paper, a new data driven dispersive analysis 
of 77 —>■ 37 / was performed. The hadronic final state 
interactions were incorporated using the Khuri-Treiman 
equation, which was solved using Pasquier inversion tech¬ 
nique. To the best of our knowledge it is the first time 
such an approach has been used in analysis of the t] de¬ 
cays. In an earlier study [3S], we illustrated the pros 
and cons of the Pasquier technique using a toy model 
with known exact solutions. The main limitation of this 
method is related to the treatment of the left-hand cuts, 
which in general are not known. We approximated them 
by a constant which is absorbed in the subtraction con¬ 
stants. As it was shown in [36], this approximation works 
very well, when the physical region does not depend 
strongly on the accurate form of the left-hand cut. On 
the other hand the advantage of the Pasquier inversion is 
that it eliminates the need for specifying the high-energy 
behavior of the absorptive parts in the physical region. 

In the analysis of the 77 —>■ Stt decays presented here, we 
have shown that with a single real parameter (gn) and 
the physical tttt partial-wave amplitudes |7] it is possible 
to reproduce the Dalitz distribution of the charged rj de¬ 
cay mode [H] . We have also verified that including more 
partial waves leads fits with comparable x^/d.o./. The 
resulting Dalitz parameters, averaged over the various 
combinations of partial waves considered in this paper 
are. 


a = 1.116 ± 0.032, 6 = 0.188 ±0.012, 
d = 0.063 ±0.004, / = 0.091 ±0.003, (41) 

g = 0.042 ±0.009. 

These are consistent, within la with the analysis of 
WASA-at-COSY having central values shifted towards 
values obtained from analysis by the KLOE Collabora¬ 
tion, which were not include in our fits. Based on the 
analysis of the charged decay we made a prediction for 
the slope parameter of the Dalitz distribution in the neu¬ 
tral decay channel, 

a =-0.022 ±0.004. (42) 

This value is above the PDG value of = —0.0317 ± 
0.0016. We speculate that the discrepancy may be a 
consequence of the WASA-at-COSY b parameter being 
significantly larger than in the earlier KLOE analysis [16] . 
We expect that in the future this issue will be resolved 
once the new KLOE data [SS] become available allowing 
a simultaneous fit of both data sets. 
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Another useful test of the amplitudes is provided by 
the ratio of neutral and charged decay rates. In the 
isospin limit this ratio does not depend on the normaliza¬ 
tion, and if the small electromagnetic isospin breaking is 
also ignored m, it depends only on the integrated Dalitz 
plot distributions. From our amplitude we find 


r (?7 —>• 37r°) 
r(?7 —>■ 7r+7r“7r*^) 


1.52 ±0.09, 


(43) 


which is consistent with the experimental value of = 
1.43 ± 0.02 [T]. We have also compared our amplitudes 
with the NLO yPT results and found the Q-value of 


Q = 21.4 ±0.4. (44) 


The error is of the statistical origin. It was computed 
through standard error propagation of the uncertainties 
arising from the tttt phase shifts, the ±3 coefficient, the 
experimental decay width r (77 —> tt+tt^tt*^) and the sta¬ 
tistical error in fitting the Dalitz plot. Inelasticity and 
higher partial waves are also potential sources of uncer¬ 
tainties EO]. 

Using the extracted Q-value and recent averages from 
the Nf = 2 ± 1 lattice computation for m = 3.42 ± 0.09 
and rUs = 93.8 ± 0.24, m we estimate the up and down 
quark masses to be 


ruu = 2.02 ± 0.14 MeV, 

rud = 4.82 ± 0.08 MeV. (45) 

The method for amplitude construction presented in 
this work can be directly applied to decays of heavier 
meson, e.g. rj' and used, for example, to test reliability of 
the isobar model. It can also be extended to incorporate 
couple-channels, which might be more relevant in decays 
of heavier mesons. 

All the material, including data and code are available 
in an interactive form online m- 


that the determined _ Qy [53 0.017 ± 0.004 is 

considerably lower than WASA-at-COSY result. It con¬ 
firms the expected correlation to the slope parameter 
in the neutral decay channel, which turned out to be 
aBESiii = _o.055 ± 0.014 ± 0.004. 


Appendix A: Isospin algebra 

In Eq. Q the isospin factors are given by 

.p(O) _ f r r 

' a /3777 ~ 2 °a/3 1 

'^a /3777 ~ 2 ~ ’ (^ 1 ) 

which satisfy. 


-pU'l Stt' 

/ j ' rjja'lS' ' ^ ’ 


(T) 


E = Kplg' , (A2) 

7/7 

' 7a/3?7' 777a^/3' ' a/3a'/3d^s^^J3^3^' * 


Here a, j5, 7, 7 / are the Cartesian isovector indices and 
isospin crossing matrices Cst and Csu, are given by 


1/3 1 5/3 


/1/3 -I 5/3' 


Cst = 1/3 1/2 -5/6 , = -1/3 1/2 5/6 

Vl/3 -1/2 1/6 / Vl/3 1/2 1/6, 


Appendix B: Kernel fnnctions 
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After submission of our manuscript a new 7 —>■ 37 r 
analysis by the BESIII Gollaboration became available 
[5]. The values for the Dalitz plot parameters (except 
the parameter /) of the charged 7 decay are compat¬ 
ible with our results within the error bars. We note 


The kernel functions in Eq. (20) are determined as 
V/Lj'L' (s, t) =2 {2L' ± 1) 


X {9{t)AjL,I'L'{s,t) — 0{—t)'SjL,I'L'{s,t)), (Bl) 


with 


POO 

^IL,I'L'{s,t) = / {C) — 

Js^(t) s' 


ds' Pis') {s'/A-ml)^ 

s-{t) 

FpL'{t)K^'{t) 

(t/4 — rn^)^' 

, ds' p{s') {s'/A — 


A 7 ds' Pis') {s'/A-mi)- 

(B2) 

u/^-'"Gr' 

and 


sUt) s- -sJ-jl{s')K-+^{s')/s' 
J^rL'{t)K-'{t) 


{t/A — mi)- 


7dPL{zs')PL'{zt), (B3) 
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FIG. 6: Integration contour C' in the complex s plane after 
Pasquier inversion. The black wiggle lines represent cuts 
attached to two branch points; {nir, ± in s-plane. The 

points labeled by a — i correspond to (a) s_(0) = —oo, 

2 2 2 2 
/, V / . Ox —ra^ , . , , o 

(b)s-(4m^) = —L-^ (c) S-(—^-2-) = 4m2, 

(d) s±((m^ - m^)^) = + m^), (e) 

2 _ 2 

s+(m,r(m^ + m^)) = (m^ - m^)^, (f) s+(4m^) = , 

(g) s+(0) = oo, (h) s+(m^(m^ - m,,)) = (m^ + m^)^, and 
(i) s+(—oo) = oo, respectively. 


where the contour C is shown in Fig. (see [3S] for more 
details). These kernel functions can be computed analyt¬ 
ically, what significantly speeds up numerical computa¬ 
tions. In the calculations presented in this paper, only 
the functions AjL jiLi{s,t) are needed and their analyti¬ 
cal representations are below in terms of 


, /-s+iO . 

(f) he') - 

Js.(t) s' 
(s' — Aml)^ 


ds' 1 


(i) s--sU(s') 


K^(s' 


PL{Zs')PL'izt), (B4) 


so that for L = 0, 


Aioj>L'(s,t) = 4 


L'-L J^I’L'(t) 


(t — 4m^)^' 


(B5) 


Ao.L'(s,t) s^h-sT 


J^io(s) 
and otherwise (L ^ 0), 


SP 


^O.L'ishKt) 


A/L,/'L'(s,t) = 4^ ^ ^L,L'(s,t) . (B6) 

The square root function U(z) is given by 

U{z) = ^(Z- (to^ - 771^)2)(z - (TO^-bTO^)2) (B7) 

in the complex z plane. Here and in what fol¬ 
lows, the phase convention for U(z) is chosen by 

C/(s ± fO) = (=F, b i) |C^(s)l for s G ((—00, (to,, — 

[(to,, — TO„)^, (to,, -I- TO,,)^], [(to,, -I- TO,r)^, oo)) respec¬ 
tively. The kinematic factor K (s)/(s p(s)) is given by the 
value of U(s) right below the two cuts attached to branch 
points s = (rriri ± to-.^.)^, i.e. K(s)/(s p(s)) = U(s — fO). 
For real s and t the physical values of Aij,/(s,t) 
correspond to the limit s -|- fO and t + iO. 


(L,LO = (0,0): 


Ao,o(s,i) — 


1 


u{s) L 


In 


Ris,t) + U(s)U(t) 


R(s,t) -U(s)U{t) 


-iTr0{ip(s,t)) 


R(s,t) = -mt + (s - TO^)(t - TO^) + TO^ (s + <), 


2 \2 


<p(s, t) = st(m„+ 3to„ -s-t)-m^(m„- to„) 


. (L,L') = (0,1): 

Ao,i(s,t) = 2 t Aa(t) + t ( 2 s + t - - 3ml) Ao.o(s,t), 


r^+it) We' 

(O') 

s_(t) 

(L,L') = (1,0): 


ds' 

Aa(i)= /(C') ^ = -ln 

Js_(t) U(s') 


(ml - ml + t + U(t)y 

4to2 t 


Ai,o(s,i) — - {Aifi(s,t) — Ai_o(0,t)) , 

Ai,o(s,t) = Ao,o(s,i) + 2 (t + to„(to,, - to,,)) A[+^(s,t) 

-b (to,, - m^)'^Ah\s,t) 

+ 2(t + mTr(mn - TO,r))(TO,, - TO,r)^Ac(s, t). 


A 


(±) 


(s,t)= he') - 
J s-(t) S 


l+h ds' 1 


1 


s_(t) S'- s H(s') s'- (to,, ± TO,,)^ 
1 


A,(s,t)= he') - 

Js.(t) s' 


s — (mjj ± TO,r)2 
l+h ds' I 


(Ao,o(s,t) - 


s_(t) s'-sU^(s') 

Ao.o(s,t) , 1 A^"^(t) 


U^(s) Am^^ mri s — (mri — m^r)^ 

1 Ah\t) 


4 TO,r ITT-rj S — (to,, -|- TO.n.)2 ’ 

r«+(b We' 1 


,,, /•«+(*) We' 

Uc') , 

Js.(t) s'-(m^z 


_(t) S - yuir^ ± TO,r)2 U(s') 

Ujt) 

TO,, (m^ ± mT^)(t ± TO,r(TO,, =F TO.,,)) 

(L,L') = (1,1): 


Ai4(s,t) — - (Ai4(s,f) — Ai^i(0,t)), 

Ai,i(s,t) = 2 Aa(t) -b4(t-bTO„(TO,, - TO„))A|^+^(t) 
-b 2 (to,, - TO„)^A^”^(t) 

-b 4 (f -b TO„(to,, — TO„))(to,, — TO„)^Ae(t) 
-b (2s -b t - TO^ - 3 to^) Ai,o(s, t), 

^(C') ^ = Ai+)(t)-A|,-)(t) ^ 

S_(t) C^^(s') 4 to „ to ,, 


/■«+(*) 
Ae(i)= he') 
J s_(t) 
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. (L,L') = (0,2): 

^o, 2 {s,t) = Af{s,t) + (2s + i - 3m^) Aa{t) 

{2s + t - ml - 3ml)‘^ - U^{t) 

+ - ^^ - - - —Ao,o{s,t), 

ds' 

= U{t) + (m^ +ml- s) Aa{t). 

. (L,L') = (2,0): 

A2,o(s, 0 = Ai(s,i) + ^ A^~\s,t) 


A^-\s,t) - A‘}-\t) 
+ 6 {t + m^ {mn - m^))— - - -^ 

3 (i + mTr{mn - TO,r))^ A^"*"^ (s, t) - a|'''^ (t) 


2 m^rriTr 


s — {mrj — mT^y 


3 (i + m^(m^ - m^))^ A^- \s,t) - A] \t) 


2 m^TOTT 

«+(*) 


s - (m^ + 


A / N njv/V s '-4 m: 

’"is? •'(»'-») c/’M 




/•«+(*) Wo' 1 1 

+\ = nr') _=_ 

^ ^ Js%) s' - s U{s') (s' - (m, ± m^) 2)2 

Ao,o{s,t) - A^^\t) Af^(t) 

(s — {mri ± m^yy s — (m^ ± m^)^ ’ 

= f (O') _ ^ _ 

JsUt) U{s'){s'-{m,±m^yy 

^4 m7r(m^ ± m^)(t ± m^(m,, =F WTr)) 

3 ( 4 m,;m ^)2 

^ U{t) / m^l/^(t) + 3 m^t(t- 4 m^) 

(p((m^ ± m^)2,i) I V3((m^ ± m^)2,t) 


/•«+(*) 

AW(i)= {C') 

J S— (t) 


nn 1 

S_(t)^ [/3(s') s'- (m^ ±7Jl;r)^ 

Ai^^(t)-Ae(t) 


4m^m7n- 

. (L,L') = (2,1): 

A2.i(s, t) = t{ 2 s + t-m ‘^-3 ml)A2fi{s, t) 

-tAe{t) + 4 :mltAc{ 0 ,t) 

+ 3 tA[. \t) + 6 t{t + m 7 r(m^ — m^)) A| ^ (t) 

01 {t -\- tWti- {fyi'q Ajyi{ty 


A„(t)= hC) 

J s_(t) 


;+(/) ds' A|+)(0-A|-)(t) 


c/5 (s') 


4 m,; mTT 


(L,L') = (1,2): 


Ai,2(s,C) — - (Ai^2(s5/) — Ai^2(0,/)) 


Ai,2(s, t) — 1 + 


2 {t + mTr{mjj — mTr)) (m^,, — m„)^ 

s — {m^i + m„)2 s — {m^ — m^^y 


2{t + m„(m,; — m^)){mjj — m,,)^ 


Ao,2(s, t) 


C/2(s) 

1 2 (/+ m,r(m,; — m„))(m,; + m„)^ 


AmT^mri 


s — {m^ + m-j^y 


1 2 {t — mTr{mrf + mTr)){mrj — m^rY 


A(+)(/) 

Al-)(/), 


AmT^mri s — {mrj — mr^y ® ’ 

/•«+(*) Ws' 1 

AW(C)= {C') —---^ — 

Js_(t) s'- (to,, ± m„)^ C/(s') 

3 /^( 2 s' + i — TO^ — 3 m'^y — U'^{t) 

X -- 

2 

= QyA^,^\t) + 6/^(to^ ± 4m,;TO„ - m^ + t)Aa{t) 
^ 3 (to 2 ± 4 m,,TO„ -ml + ty -U^ {t) ^ 

+ 2 
S_L\ /•«+(*) Wo' 

_L ™ ^2^ 


,,, /-s+l.'; Wo' 

a! ^(C) = / (C') 77;^ (s' - ("ir/ ± 

a/ S_ (t) ^ ) 

= C/(/) =F 2 TO„TO,, Aa(/). 


(L,L') = (2,2): 


A ' ,,_3C2(2s + t-TO2-3TO2)2-C/2(t)^ ^ 

A2,2(s,C) — -^-A2^o(s,/) 

+ 6/^A„(s,t) + 6i^(2s + / - TO^ - 3ml)Ao{t), 
/•s+d) 1 si 

A ^ n ^ ^ 




3s'^{2t + s' — TO^ — 3m’yy — U‘^{s') 

^ 2 

= ^ A|,"^(C)-(s-(TO,;-TO„)^)A^"^(t) 

+ 6 (/ + mTr{m,j — TO„)) 

X Al-'W-(S-K + .„W^)^1<|X^^ 

„ >a[aW(()-A.(() 


+ 6 (/ + mrr{mn — to„)) 


4 m,, mrr 


+ {s- (to,, - mrrf) 


2 A,(t)-A^,y\t)-Ai-\t) 

(4 m,, m^n-)^ 


Ai^^(C) - (s - (to,, - to„)^)A, 


+ 2 m^ [Ae(/) - sAc( 0 ,t)] 
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/■«+(*) 1 1 

A,(s,t)= J{C') 


3s^ {2t + s'— — 3m^y — U^{s') 


= l^fc \t)+Q{t + m^{mr, - 

2 Amr^niT:- 

^,2AS)-Ai^\t)-Ai-\t) 


— 6 {t + m7r(TO^ — TO^)) 

- \ (AeW -4m2 Ac(0,i)) . 


(4m^ m^) 
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